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ABSTRACT 

We reinvestigate the correlation between black hole mass and bulge concentration. With an 
increased galaxy sample (totalling 27), updated estimates of galaxy distances, black hole masses, 
and Sersic indices n — a measure of concentration — we perform a least-squares regression 
analysis to obtain a relation suitable for the purpose of predicting black hole masses in other 
galaxies. In addition to the linear relation, log(Afbh) = 7.81(±0.08) + 2.69(±0.28) log(n/3) with 
Eintrinsic = O.SI^q^q^ dcx, wc investigated the possibility of a higher order M^i^-n relation, finding 
the second order term in the best-fitting quadratic relation to be inconsistent with a value of 
zero at greater than the 99.99% confidence level. The optimal relation is given by log(Afbh) = 
7.98(±0.09) + 3.70(±0.46)[log(n/3)] - 3.10(±0.84)[log(n/3)]2, with eintrinsic = 0.18to:o6 dex and 
a total absolute scatter of 0.31 dex. Extrapolating the quadratic relation, it predicts black holes 
with masses of ~ lO^M© in n = 0.5 dwarf elliptical galaxies, compared to ~ IO^Mq from the 
linear relation, and an upper bound on the largest black hole masses in the local universe, equal 
to l.2tli X IO^Mq. 

In addition, we show that the nuclear star clusters at the centers of low-luminosity elliptical 
galaxies follow an extrapolation of the same quadratic relation — strengthening suggestions for a 
possible evolutionary link between supermassive black holes and nuclear star clusters. Moreover, 
we speculate that the merger of two such nucleated galaxies, accompanied by the merger and 
runaway collision of their central star clusters, may result in the late-time formation of some 
supermassive black holes. 

Finally, we predict the existence of, and provide equations for, an Afbh-^o relation, in which 
fiQ is the (extrapolated) central surface brightness of a bulge. 

Subject headings: black hole physics — galaxies: bulges — galaxies: fundamental parameters — galaxies: 
structure 



1. Introduction 

While perhaps not surprising to AGN-astronomers, 
over the past five to ten years, the notion that a 
supermassive black hole (SMBH; M^h > IO^Mq) 
resides at the heart of every significantly large 
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(Mbuige > IO^Mq) galaxy bulge, even inactive 
bulges, has changed from a dubious idea to a 
mainstream belief. Support for the tide-of-opinion 
change has, in part, arisen from studies of our own 
galaxy. Schodel (2003) and Ghez et al. (2005, and 
references therein, see also Broderick & Narayan 
2006) have shown that the mass required inside of 
the innermost resolved volume of the Milky Way 
is sufficiently large that it rules out alternatives to 
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a SMBH, such as a cluster of dead stars or stellar 
mass black holes which, if they did once exist, have 
surely now merged to form a single, massive ob- 
ject (e.g., Miller 2006). Although there is only one 
other sufficiently well resolved galaxy where such 
a conclusion can be drawn, NGC 4258 (Miyoshi 
et al. 1995), it seems reasonable to accept that 
the dark concentrations of mass at the centers of 
other galaxies are also SMBHs, and we adopt this 
convention, or at least terminology, here. 

Roughly some three dozen galaxies are close 
enough that their central, SMBH mass has been 
measured directly through its influence on the mo- 
tion of the surrounding gas and stars (Kormendy 
& Gebhardt 2001; Merritt fc Ferrarese 2001a; Fer- 
rarese & Ford 2005, their Table II) ^ To acquire the 
mass of the central dark object in other (inactive) 
galaxies requires an alternative approach, and a 
number of indirect means to do so have been pro- 
posed (e.g., Ferrarese & Merritt 2000; Gebhardt 
et al. 2000; Marconi & Hunt 2003; Novak, Faber 
& Dekel 2006). 

Perhaps the most often used relation — due to 
the small level of scatter (~0.3 dex) — is the Mbh- 
(7 relation (Ferrarese & Merritt 2000; Gebhardt 
et al. 2000). With an equivalent level of scatter, 
and the advantage of requiring only galaxy images 
rather than spectra, is the M^h-n relation (Gra- 
ham et al. 2001, 2003a), where n is a measure of 
the concentration of the stars within the bulge. 
More precisely, n is the inverse exponent from 
the best-fitting Sersic (1963) i?^/" light-profile (see 
Graham & Driver 2005 for a review of this model). 
Other relations, often reported to have more scat- 
ter than the above two relations, are the Afbh L 
and the Mbh Mbuigc relations (Kormendy 1993; 
Magorrian et al. 1998). Although, recent stud- 
ies which have excluded the disk-dominated spiral 
galaxies (e.g., McLure & Dunlop 2002), or used 
only the bulge luminosity from the disk galaxies 
after performing an accurate bulge/disk decom- 
position (e.g., Erwin et al. 2002; Marconi & Hunt 
2003) have obtained a relationship with a similarly 
low level of scatter (see also Haring & Rix 2004) . 
All of these empirical relations can be, and indeed 
must be, used to constrain any complete theory 

We do not use the eight galaxies listed in part two of Ta- 
ble II from Ferrarese & Ford 2005 for which the SMBH 
masses might be in error. 



or model of galaxy/SMBH co-evolution. Further- 
more, all of these relations can be used to gauge 
the mass of SMBHs in other galaxies. 

In this paper we use the Sersic index n, together 
with updated SMBH masses (Section 2), to de- 
rive the first Afbh n relation constructed for the 
purpose of predicting Mbh in other galaxies (Sec- 
tion 3). This relation may differ from the intrinsic, 
astrophysical relation presented in Graham et al. 
(2003a) due to the different method of regression 
that is required. Armed with such a relation, one 
requires only images — which need not even be 
photometrically calibrated — to predict accurate 
SMBHs in other galaxies. Additional advantages 
with the use of a global measure such as n, are 
that it is not heavily affected by possible kine- 
matical sub-structure at the center of a bulge, nor 
by rotational velocity or vertical dispersion of an 
underlying disk, nor aperture corections. Further- 
more, the quantity n is cheap to acquire in terms 
of telescope time and unlike absolute magnitudes 
and masses, it does not depend on galaxy distance 
nor an uncertain mass-to-light ratio. 

As noted in Graham et al. (2001), there is no a 
priori reason to presume that the relation between 
log Mbh and logn is linear, and we therefore ex- 
plore the suitability of a a quadratic equation. In 
so doing, we find the second order term is inconsis- 
tent with a value of zero at the 99.99% confidence 
level (Section 3.2). Implications at the low- and 
high-n end of the relation are discussed, relative 
to the optimal linear fit. 

In a forthcoming paper we will apply this 
quadratic relation to the Millennium Galaxy 
Catalog (Liskc ct al. 2003; Cross ct al. 2004; 
Driver et al. 2006) containing 10 095 galaxies 
— which have been modeled with a Sersic-bulge 
plus exponential-disk (Allen ct al. 2006) to de- 
termine the local supermassive black hole mass 
function (e.g., Salucci et al. 1999; Yu & Tremaine 
2002; Granato ct al. 2004; Shankar ct al. 2004) 
and space density in both early- and late-type 
galaxies. 

In Section 4 we expand the Mb^-n diagram into 
an M^cQ-n diagram, where Mmco is the mass of 
the central compact object, which may be a SMBH 
or a nuclear star cluster. We show that the nuclear 
star clusters in early-type galaxies appear to fol- 
low the curved Mbh^ relation defined by galaxies 
with SMBHs, and we discuss some of the implica- 
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tions this may entail. In Section 5 we show that 
our Mbh-?^ relation is independent of the Hubble 
constant. 

Given there are now several reflations between 
SMBH mass and the properties the host bulge, in 
Section 6 we briefly present some musings as to 
what may be the fundamental primary relation. 
In this section we derive a new set of equations re- 
lating the SMBH mass to the (extrapolated) cen- 
tral surface density of the host bulge. Finally, in 
Section 7 we provide a brief summary of the paper. 

2. Data for the M^h—n relation 

Our sample is comprised of the 21 galaxies used 
in Graham et al. (2001) plus an additional six new 
galaxies. 

A discussion of the original galaxy light-profiles 
can be found in Erwin ct al. (2002), sec also Tru- 
jillo et al. (2004). Due to updated galaxy dis- 
tances and SMBH mass measurements, and a re- 
fined analysis of the major-axis light-profiles, some 
of the 27 data points (14 E and 13 SO/Sp) given 
in Table 1 are slightly different to those published 
in Graham et al. (2001, their table 1) and shown 
in Graham et al. (2003a, their figure 1). However, 
only for two galaxies has the Sersic index changed 
by more than 20% the typical uncertainty on 
this index (e.g., Caon, Capaccioli, & D'Onofrio 
1993). We briefly comment on these two galax- 
ies here and present their new Sersic flts in Ap- 
pendix A. 

Although the bulge of the Milky Way is nowa- 
days recognised as having an exponential (n = 1) 
light-profile, we have fitted the near-infrared data 
from Kent, Dame, & Fazio (1991) and report an 
index of n = 1.32 rather than exactly 1. The 
second galaxy in question is NGC 4564. Sub- 
sequent to the analysis in Trujillo et al. (2004), 
NGC 4564 is now recognized as an SO galaxy that 
was mis-classified as an E galaxy. An iii/"-bulge 
plus exponential-disk decomposition has therefore 
been performed. This resulted in its bulge Sersic 
index n increasing from 2.1 to 3.2. 

Four of the six new galaxies are listed in Fer- 
rarese & Ford (2005) (NGC 3115, 4486, 4649, 
4697) and our Sersic fits to their light-profiles are 
also presented in Appendix A. In addition, we have 
used the bulge/disk decomposition shown in Gra- 
ham (2002) for the supposedly "compact Ellipti- 



cal" galaxy NGC 221 (M32), which was not in- 
cluded in our original paper. We have also been 
able to include NGC 1399 (Houghton et al. 2006), 
which previously had no direct SMBH mass esti- 
mate. This galaxy's light-profile has been modeled 
in D'Onofrio, Capaccioli, & Caon (1994) and we 
adopt their Sersic index having obtained the same 
value from our own fitting of their data. 

Despite our attempt to include all of the low- 
est and highest mass SMBHs, as these data points 
can have substantial weight on any fitted relation, 
we have not included the peculiar elliptical galaxy 
IC 1459 due to uncertainty on both its SMBH 
mass and its Sersic index n. This galaxy displays 
clear signs of past interaction, evidenced by stellar 
tidal tails (Malin 1985) and stellar shells and rip- 
ples at large radii (Forbes & Reitzel 1995). Due 
to its unrelaxed, disturbed morphology, its light- 
profile is not well fitted with an R^^" model. It 
also possesses a fast counter-rotating stellar core. 
While the stellar dynamics of the core suggest a 
SMBH mass of 2.6 x lO^M©, the gas dynamics 
reveal the mass could be as low as 3.5 x lO^M© 
(Cappellari et al. 2002). 

The upper section of Table II in Ferrarese 
& Ford (2005) Hsts 25 galaxies with SMBH 
masses derived from resolved dynamical stud- 
ies. A further five galaxies in this list have 
marginally resolved mass estimates, with 0.39 < 

^sphere of infiucncc/^rcsolvcd ^ 0.92. Aside frOm 

IC 1459, we are only missing quality i?-band im- 
ages for three of these 30 galaxies: NGC 3608, 
NGC 5128 (Centaurus A) and the relatively dis- 
tant {z = 0.056) galaxy Cygnus A. 

Most of our 27 galaxies have distances from 
Tonry et al. (2001), and updated SMBH mass es- 
timates have come from Tremaine et al. (2002), 
with exceptions noted in Table 1. Three galax- 
ies in our sample do not have 'Surface Bright- 
ness Fluctuation' distance measurements in Tonry 
et al. (2001). For the relatively distant galaxies 
NGC 6251 and NGC 7052, we used their heliocen- 
tric velocities (Wegner et al. 2003) and a Hubble 
constant Hq = 73 km s"-'^ Mpc^^ (Blakeslec et 
al. 2002) — consistent with the HST Key project 
value of 72±3 ± 7 (Frcedman et al. 2001) and the 
WMAP value of 73^1 (Spcrgel et al. 2006). For 
the remaining Virgo cluster galaxy NGC 4342, we 
used the mean Virgo cluster distance of 17.0 Mpc 
(Jerjen, Binggeli, & Barazza 2004). The above 
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three galaxies have shghtly different SMBH masses 
in Table 1 to those in Tremaine et al. (2002) due 
to the slightly different distances used. 



Table 1 
Galaxy Parameters 



(jrcilaxy 


Distance 
Mpc 


(lO^Mo) 


n 




Elliptical Galaxies 




NGC 821 


24.1 




4.00tS1? 


NGC 1399'^ 


20.0 




16.8+2-8n 


NGC 3377 


11.2 


1.00+0-? 


3.04+0-6} 


NGC 3379 


10.6 


-, oc+0.73 
^■'^^-0.73 


4.29+°-?^ 


NGC 4261 


31.6 


5.20+}-; 


7.30+}-^! 


NGC 4291 


26.2 


3.10+°i 


4.02to;«° 


NGC 4374 


18.4 


4 64+3-46 




NGC 4473 


15.7 


1 in+°-'^° 

-'-•-1-U_Q 79 


' "5 n /ifi 


NGC 4486 


16.1 


34.3+9-7 


6.86l}-37 


NGC 4649 


16.8 


2o.ot^:;5 


6.04+}-2i 

— i .uu 


NGC 4697 


11.7 


1.70+0-2 


4.00+0-80 

— U.u ( 


NGC 5845 


25.9 


2.40t?-4 


o 99+0.64 
54 


NGC 6251 




5.80+2-0 


11.8+2-36 


NGC 7052 




3.40^2;4 


4.551°;?^ 




Bulges of Disk Galaxies 




Milky Way*" 


0.008 




1 00+0.26 
l.oz_o 22 


NGC 221^= 


0.81 


025+0-005 


-1 r;-|+0.30 


NGC 1023 


11.4 


44+0-0-5 


2 m+0-40 
^•'Jl-0.34 


NGC 2778'* 


22.9 


14+0-08 


1 60+0-32 
l.DU_o.27 


NGC 2787^ 


7.5 


0.411°:°^ 


-0.33 


NGC 3031 


3.9 


0.68^n^ 


Q ofi+0.65 
"^•^0-0.54 


NGC 3115-'' 


9.7 


9 20+3-00 


13.0+2-60 


NGC 3245 


20.9 


2 in+0-05 


4 -31+0.86 
+ "Jl-0.72 


NGC 3384^ 


11.6 


16+0-01 


1 79+0.34 
^^•'^-0.29 


NGC 4258S 


7.2 


OQ + 0.01 


2.04t°:^i 


NGC 4342^ 


17.0 


3.30+}-? 




NGC 4564'^ 


15.0 


^fi+O-03 


q 1 r+0.63 


NGC 7457 


13.2 


n 035+0-011 
u.uoo_o 014 


1 oq+0.37 
-'-•°'J-0.31 



"^B-band image. 



''2.4-^m minor-axis profile from Kent, Dame, & Fazio (1991). 
'^Taken from Graham (2002). 

<^NGC 2778 is a misclassified SO galaxy (Rix, Carollo, & Freeman 1999). 
<^HST F814W image. 
^7-band image. 



5 



^Thuan-Gunn r image. 

'^NGC 4564 is a misclassified SO galaxy (Trujillo et al. 2004). 

Note. — Distances are taken from Tonry et al. (2001, their table 1), except for the Milky Way (Eisenhauer et al. 2005), 
NGC 4342 (Virgo cluster distance modulus = 31.15, Jerjen, Binggcli, & Barazza 2004), NGC 6251 (wcmb=7382 km s'^, 
Wcgncr ct al. 2003), and NGC 7052 (?;cmb=4411 km s'^, Wegner et al. 2003). These four galaxies are not listed in Tonry 
et al. (2001). A Hubble constant of Hq = 73 km s'^ Mpc'^ (Blakeslee et al. 2002; Spergel et al. 2006) has been used for 
the latter two galaxies. The SMBH masses are from the compilation in Tremaine et al. (2002), except for the Milky Way 
(Ghez et al. 2003, see also Beloborodov et al. 2006), NGC 821 (Richstone ct al. 2006), NGC 3379 (Gebhardt et al. 2000; 
see also Shapiro et al. 2006), NGC 4486 (Macchetto et al. 1997) and NGC 3115 (Emsellem, Dejonghe, & Bacon 1999). 
Our sample includes a further three galaxies not listed in Tremaine et al. (2002). The SMBH mass for NGC 3031 and 
NGC 1399 are from Merritt & Fcrrarcse (2001a) and Houghton ct al. (2006), respectively, and the mass for NGC 4374 is 
from Maciejewski & Binncy (2001, with updated errors taken from Kormendy & Gebhardt 2001). The Sersic indices are 
from the major-axis, light-profiles of ii-band images, except where noted otherwise. 
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3. The Mbh— ^^ correlation 
3.1. A linear relation 

The Mbh-n relation presented in Graham et 

al. (2003a) was constructed using a bisector Hnear 
regression analysis which treated both variables 
equally. In that study, the relation was derived 
with the goal of determining the intrinsic physical 
relation between Mbh and the Sersic shape param- 
eter n of the bulge. Here our objective is different 
because we wish to obtain a relation that can be 
used to predict values of Mbh in other galaxies from 
their observed bulge Sersic index n. We there- 
fore desire an M^^i^-n relation which minimizes the 
scatter in the quantity to be predicted, and have 
thus performed an ordinary least squares (OLS) 
regression of Mbh on n for those local galaxies for 
which both quantities are known (Table 1)^. 

Following Graham et al. (2003a), we have used 
a 20 per cent measurement error on the values of 
n, or more specifically, we have assigned an error 
of ±log(1.2) to log(n). Later on we explore the 
influence of varying this quantity. Several factors 
can contribute to the size of this term, including 
errors in the sky-subtraction, uncertainties in the 
point-spread-function, the presence of bars which 
are typically not modeled in bulge/disk decompo- 
sitions, the influence of additional nuclear compo- 
nents such as star clusters or nuclear disks, and 
partially depleted cores. The latter two issues can 
be dealt with by either simultaneously fitting a 
Sersic function plus some additional function to 
account for the excess flux above that of the host 
galaxy (e.g., Graham & Guzman 2003; Ferrarese 
et al. 2006a), or with the use of the core-Sersic 
model (Graham et al. 2003b; Trujillo et al. 2004). 
However, due to these various issues, it can be dif- 
ficult to acquire reliable individual uncertainties 
on the Sersic index for every galaxy, hence our use 
of an average relative error. Similarly, the use of 
a fixed relative error of 13% and 5% was assigned 
to the velocity dispersion term used by Merritt 
& Ferrarese (2001b) and Tremaine et al. (2002), 
respectively, in their construction of the Mbh-c 
relation (with the exception that they both used 
a 20% uncertainty for the Milky Way's velocity 
dispersion). 

^Scc Fcigclson & Babu (1992) for a clear exposition behind 
the rationale of when to use what type of regression. 



We have used Tremaine et al.'s (2002) modified 
version of the routine FITEXY (Press et al. 1992), 
to solve the equation y = a + bx, by minimising 
the quantity 

,.2 _ {yi-a- bxjf 

USy^+bHx^ + e^- ^' 

The measurement errors on Xi and yi are de- 
noted by Sxi and Syi, and the intrinsic scatter e is 
searched for by repeating the fit until /{N — 2) 
equals 1. The uncertainty on e is obtained when 
the reduced chi-squared value, / {N — 2), equals 
1 ± ^J2/N. Doing so, one obtains 

log(Mbh) = 2.69(±0.28)log(n/3) + 7.81(±0.08), 

(2) 

with e = 0.3ll^:^? dex in log Mbh- This fit is 

shown in Figure 1. The total absolute scatter in 
log Mbh is 0.39 dex. Using either an optimistic 
~10% estimate for the uncertainty on n (specif- 
ically, using log[n] ± log[l.l]), or a more liberal 
~25% uncertainty (we used \og[n\ ±log[1.25]), see 
for example Caon, Capaccioli, & D'Onofrio (1993), 
had no significant (l-(j) affect on cither the slope 
or intercept of the above relation. The new intrin- 
sic scatter was 0.35 dex and 0.27 dex, respectively. 

The maximum l-rr cirror on the predicted value 
of Mbh in galaxies for which n is known is acquired 
by assuming uncorrected errors on n and the slope 
and intercept of the Mbh-w relation. Gaussian 
error propagation for the linear equation y = (b± 
5b) [x ± 5x) + {a± 5a) , gives an error on y equal to 

5y = 

^/{dy/db)'^(5b)-^ + idy/da)'^i5a)'^ + {dy / dx)"^ {5x)'^ 
= \/a;^((5&)2 + (5a)2 + &2(5a;)2_ 

In the presence of intrinsic variance, e, the uncer- 
tainty on y will be greater, such that, assuming 
the intrinsic variance is in the y coordinate, 

5y = ^x^{5b)^ + (<5a)2 + b'^{5x)^ + e^. 

For our expression (equation 2) we have x = 
log(n/3), and so dx/dn = l/[ln(10).n], and there- 
fore 

(51ogMbh)' = [log(n/3)]2(0.28)2 + (0.08)2 
+ [2.69/ln(10)]2[5n/n]2 + (0.31)^. (3) 
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Fig. 1. — Correlation between a galaxy's supermassive black hole mass and the shape parameter (i.e. Sersic 

index n) of its dynamically hot component. The Pearson linear correlation coefficient r is given, as is the 
Spearman rank-order correlation coefficient r^. (The uncertainties on the data points were not used when 
computing these correlation coefficients.) The regression line shown in the left panel was obtained using a 
modified version (Tremaine et al. 2002) of the routine FITEXY (Press ct al. 1992, their Section 15.3), see 
equation 2. Consistent results were obtained using the ordinary least squares (log(Mbh) 1 log(n/3)) linear- 
regression routine from Akritas & Bershady (1996), see equation 4. The middle panel shows the Ax^ = 1-0 
and 2.3 boimdaries aroimd the optimal intercept, a = 7.81, and slope, b = 2.69. The projection of the 
Ax^ = 1.0 ellipse onto the vertical and horizontal axis gives the 1-a uncertainties 5a and db, respectively. 
The Ax^ = 2.3 ellipse denotes the 1-a two-dimensional confidence region. This has been mapped into the 
right panel, and is traced by the two solid curves. The dashed lines in this panel arc the (more commonly 
used) approximations obtained using (a ± 6a) and (b ± 6b). The two confidence regions agree well, although 
the region traced by the dashed lines is, as expected, smaller. 



Given the debate over the slope of the Afbh f 
relation (Merritt & Ferrarese 2001b; Tremaine et 
al. 2002; Novak et al. 2006), it is of interest to 
know how much the uncertainty in the slope and 
intercept of the M^i^-n relation may contribute to 
the uncertainty in the predicted SMBH masses. 
From equation 3, if one measures a bulge to have 
n = 3 with 6n/n = 0.2, then the error in n con- 
tributes 96% of 6 log Mbh- That is, the uncertainty 
on the slope and intercept of equation 2 are not 
substantial contributors to the error budget on 
log Mbh- If one has a galaxy with n = 1 or 9, 
and again 6n/n = 0.2, then the combined error 
from the uncertainty on the slope and intercept 
contributes only 14% of the error on log Mbh- 

We have additionally used the OLS regression 
analysis BCES(logMbh | logn) from the code of 
Akritas & Bershady (1996), which allows for both 
measurement errors and intrinsic scatter. For our 
sample of iV = 27 galaxies, we obtained 

log(Mbh) = 2.68(±0.40) log(n/3) + 7-82(±0-07), 

(4) 

with e = 0.30^0 07 dex in log Mbh- Using an uncer- 
tainty of ±log(l-25) or ±log(l.l) for the value of 



log(n) had no significant (Icr) affect on cither the 
slope or intercept of the above relation — which 
agrees well with that in equation 2. 

3.1.1. Symmetrical regression 

To obtain the intrinsic astrophysical relation, 
the modified version of FITEXY which minimizes 
equation 1 should not be used. The reason is be- 
cause it is biased — to produce a low slope — by 
the minimisation of the intrinsic variance, e, along 
the y (i.e. log Mbh) axis. If the minimisation is in- 
stead performed along the x (i.e. logn) axis, then 
the term in the denominator of equation 1 will 
be replaced with 6^e^. Making this substitution, 
and performing the new regression'^, one obtains 
b = 3.10 ± 0-33, a = 7.78 ± 0-09, and e = 0-11^^2 
dex in logn. 

The average of the above two slopes from the 
modified FITEXY routine, (3.10-|-2.69)/2 = 2.90, 
agrees well with the slope obtained using the bisec- 
tor linear regression routine BCES from Akritas & 



We have ignored unknown, but possible, selection bound- 
aries in logn that could bias such a fit (see, e.g., Lynden- 
Bell et al. 1988, their Figure 10). 
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Bershady (1996), which gives 

log(Mbh) = 2.85(±0.40)log(n/3) + 7.80(±0.07). 

^ (5) 

This relation is in good agreement with the (intrin- 
sic astrophysical relation) presented in Graham et 
al. (2003a). 

3.2. A curved relation 

As noted in Graham ct al. (2001, 2003a), we 
have no a priori knowledge that the Afbh re- 
lation is linear. For this reason, in those pa- 
pers we employed the use of both linear and 
non-linear statistics to measure the correlation 
strength. Here we go one step further by fitting a 
quadratic to the data. 

Wc stress that, from a physical stand point, we 
do not know what the actual form of the rela- 
tion should be. The quadratic equation which we 
adopt is an empirical model. In a Taylor scries 
expansion it is simply the next order term. The 
Mi,h-n data may in fact be described by a double 
power-law, however this would require the use of 
four free parameters (a low- and high-mass slope, 
and a transition mass and transition Sersic index). 
The quadratic relation has only three parameters 
and is the adopted model for explorations of non- 
linearity in the Mbh-c data (Wyithe 2006a,b). 

In passing we note that there have been c;laims 
that the Mbh-c relation may not be linear, but 
has either negative curvature (e.g., Granato et al. 
2004; Cirasuolo ct al. 2005) positive curvature (e.g. 
Hopkins & Hernquist 2006), no curvature at even 
the 0.75-1.5 a level (Wyithe 2006b)*, or curves 
down at low SMBH masses and up at high SMBH 
masses (Sazonov et al. 2005). Given that the L-a 
relation is not linear, having a slope of ~4 at the 
bright end and ^2 at the faint end (Tonry 1981; 
Davies et al. 1983; Held et al. 1992; De Rijcke et 
al. 2005; Matkovic & Guzman 2005), and «/ the 
M\y\i-L relation is linear, then the M^i^-cr relation 

^The tighter constraint of 0.75 a (that is, the factor in front 
of a quadratic term is inconsistent with zero at only the 0.75 
(T level), comes from using SMBHs with resolved sphere's 
of influence. Including galaxies with unresolved sphere's of 
influence and adding in single epoch reverberation mapping 
masses (uncertain to factors of 3-4) weakens the result of 
a purely linear Mbh-c relation, with the probability that 
the second order term does not equal zero ruled out at the 
(still weaJf) 1-5-2 a level. 



obviously cannot be linear, but must have a posi- 
tive curvature. Of course, the Mbh-i relation may 
not be linear. 

We fit the quadratic equation y = a + bx + cx^, 
to the (Mbh,»^) data by minimising the statistic 

2 ^ {yi-a-hxj-cxjf 
^ ^,6yi^ + {2cxi + b)^Sxi^ + e^' ^' 

where = log(Mbh,i), Xi = log(nj/3), and e is the 
intrinsic variance which, given our objective of pre- 
dicting new SMBH masses, we attribute entirely 
to reside in the y (log Afbh) direction. Solving for 
xV(A^-3) = 1, we find 

log(Mbh) = 7.98(±0.09)-H3.70(±0.46)log(n/3) 
-3.10(±0.84)[log(n/3)]2, (7) 

with e = O.lSl^Q Q6 ^^'^ total absolute scatter 
in logMbh is 0.31 dex. From Figure 2, the term c 
is inconsistent with a value of zero at greater than 
the 99.99% confidence level'''. 

With small data sets, formal errors can under- 
estimate the true error on a fitted parameter. 
To explore the danger that we may be under- 
estimating the uncertainty on the parameter c in 
front of the quadratic term in equation 7, we use 
Bootstrap Monte Carlo simulations. The random 
number generator rani from Press et al. (2002) 
was used in a bootstrapping process that involved 
sampling (with replacement) from the original 27 
data points. One thousand data sets containing 
27 data points each were then individually fitted 
to find the optimal log-quadratic equation, ex- 
actly as done with the original data. This then 
gave 1000 new estimates for each of the three pa- 
rameters (a, 6, c). The range in values covered by 
the central 68.3% of these three data sets gives 
one an estimate of the l-a confidence intervals 
— without making any assumption about (the 
Gaussianity of) the distribution. This parame- 
ter range is found to be (7.88.. .8.04, 3.26.. .4.24, 
—3.98...— 2.30). This is in good agreement with 
the l-a parameter uncertainties obtained from 

® Excluding the two galaxies with the highest SMBH masses, 
the total absolute scatter reduces to 0.25 dex, the intrinsic 

scatter drops to zero, and one finds that a = 7.90, h = 
3.22, and c = —2.55, only a l-a deviation from the values 
obtained using the full data set. Using the full data set, 
but with the intrinsic scatter set to zero, the value of c is 
still inconsistent with a value of zero at the 3-<t level. 
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Fig. 2. — Left panel: The data points are the same as in Figure 1, but this time the optimal 3-parameter 
log-quadratic relation (equation 7) is shown. The absolute scatter in logAfbh is 0.31 dex. The shaded 
area corresponds to the Ax^ = 3.53 (l-cr) confidence region. Middle and right panel: The inner, central 
and outer solid contours denote the Ax^ = 1.00, 9.00 and 15.1 boundaries, respectively: their projections 
onto the vertical and horizontal axes give the 68.3% (l-o") 99.73% (3-(t) and 99.99% confidence intervals 
on the individual parameters in equation 7. The inner and outer dashed lines denote the Ax^ = 2.30 
and 11.8 boundaries, respectively, delineating the 1-a and 3-a two-dimensional confidence regions. The 
dotted lines show the Ax^ = 3.53 and 14.2 boundaries, respectively. They correspond to the 1-a and 3-cr 
three-dimensional confidence ellipsoids, (see Press et al. 1992, their Section 15.6). 



Figure 2 and given in equation 7. We therefore 
conclude that we are not under-estimating our er- 
rors, and that the Afbh ^ relation is indeed curved. 

Using equation 7, and the value of n in other 
galaxies, the maximum error on the predicted 
value logMbh, denoted by 51ogMbh, can be writ- 
ten as 

{6yf = x^idcf + x^iSbf + (Sa)^ 
+ {2cx + by{6xy + e^, 

which equates to 

{S log Mbh)' = [log(n/3)]^ + [log(n/3)] V4 + (0.09)^ 
-|-[3.70 - 6.201og(n/3)]^[5n/n]V[ln(10)]^ -|- (0.18)^(8) 

Assuming 5n/n = 0.2, Table 2 gives the value of 
(5 logMbh for different values of n. Obviously the 
error in the slope of the relation (similarly for the 
linear relation) causes the masses to be less well 
constrained at the end of the relation. Due to the 
steepness of the quadratic relation at low n, the 
uncertainty on the estimated SMBH masses can 
be large there. 

In passing we note one of the implications of a 
curved M^h-n relation. Either the Mbuige-?^ rela- 
tion must be curved, or the Mbh/Afbuige ratio can 
not be constant with mass. 



3.3. The high-mass end 

An obvious consequence of our second order, 
log-quadratic, relation is the different expected 
masses at the ends of the relation. For example, 
at the high mass end, we no longer predict in- 
finitely large SMBHs. Obviously a linear M\^\i-n 
relation or a linear A/bh <^ relation, and a posi- 
tively curving log-quadratic Afbh o" relation, im- 
ply infinitely massive SMBH masses at the high-cr 
end. In contrast, our negatively curving Afbh n 
relation suggests a maximum mass limit to which 
the universe has built SMBHs^. The maximum 
mass occurs where the derivative of equation 7 
equals zero, i.e., where dlog(Afbh)/(ilog(n/3) = 
2ca;-|-6 = 2(-3.10)log(n/3)-h3.70 = 0. The value 
of n at this rather broad peak and turnover is 11.9. 
Scanning the Ax^ — 3.53 (l-cr) confidence region 
shown in Figure 2, so as to provide error bounds, 
this peak corresponds to a maximum SMBH mass 
of 1.21q;4 X lO^Af0. This upper value matches weU 
with the uppermost mass reported in several ob- 
servational studies (e.g., McLure et al. 2004; Xie, 
Zhou, & Liang 2004; Plan, Falomo, & Treves 2005; 
Sulentic et al. 2006; Kollmeier et al. 2006) and 

'We do not equate this with a upper mass to which the 
universe can build SMBHS, only with what the universe 
has built. 
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Table 2 

Black hole mass uncertainties 



;) = 


0.5 1.0 2.0 .'-S.O 


4.0 


8.0 


10.0 




Linear Fit (equation 


2) 






log Mhh 


5.72 6.53 7.34 7.81 


8.15 


8.96 


9.22 


(51ogMbh 


0.45 0.42 0.40 0.40 


0.40 


0.41 


0.42 




Quadratic Fit (equation 7) 






log Mbh 


3.22 5.51 7.23 7.98 


8.39 


8.99 


9.07 


51ogMbh 


1.05 0.70 0.47 0.38 


0.33 


0.36 


0.43 



Note. — Uncertainty, (5 log Mbh dex, on an estimated black 
hole mass using the linear and quadratic M^h-n relations as 
a function of n. The steepness of the quadratic relation at 
low n results in both a better fit to the data for 1 < n < 2, 
but also an increased uncertainty on the predicted value of 
log Mbh. 



intriguingly with the upper mass limit for axion 

bubbles (Svidzinsky 2006). We note that this peak 
is reached somewhat asymptotically as the value 
of n increases, and that we do not expect galaxies 
with larger values of n to have significantly smaller 
SMBH masses than the peak mass. 

We additionally note that some studies of high- 
rcdshift quasars have reported SMBH masses up 
to IO^^Mq. If correct, such objects would repre- 
sent a challenge for the quadratic relation defined 
by the current local sample. Of course, to properly 
solve this discrepancy requires direct SMBH mass 
measurements. With the upcoming ten-fold in- 
crease in spatial resolution (compared to the Hub- 
ble Space Telescope), the next generation of ex- 
tremely large telescopes, such as the Giant Mag- 
ellcn Telescope, will be able to resolve the sphere 
of infiuence around SMBHs in galaxies ten times 
further away. The resultant 1000-fold increase in 
survey volume should help address the issue of just 
how massive SMBHs are. 

It is also worth remarking that the upper mass 
limit we find is somewhat of a consequence of our 
use of a quadratic relation. That is, had we used 
a four-parameter double power-law, we would not 
observe such behaviour. Although, Figure 2 does 
suggest that some kind of asymptotic function. 



which reaches a maximum SMBH mass, would also 
provide a good description to the data. We will, 
however, postpone a comparison of assorted arbi- 
trary functions for a later time. 

The two galaxies with the highest SMBH 
masses (NGC 4486 and 4649) appear as outliers 
from the curved relation. We have checked their 
Scrsic indices are correct and it may be of value 
to re-examine their black hole masses. While this 
is beyond the scope of the present paper, we note 
that there is reason to suspect that NGC 4486 
(M87) might have had its black hole mass over- 
estimated. Maciejewski & Binney (2001) have 
shown that the SMBH mass in NGC 4374 was re- 
duced by a factor of nearly four, from 1.7 x lO^M© 
to 4.64x 1O^M0, after the slit- width was appropri- 
ately deah with. NGC 4486, like NGC 4374, has 
similarly had its SMBH mass derived from the 
emissions of a circumnuclear gaseous disk using 
slit-spectroscopy. These galaxies also appear to 
have overly large SMBH masses when compared 
with the theoretical models of Menci et al. (2006, 
their Figure 1). It may be advantageous to con- 
firm this galaxy's SMBH mass using AO-assisted 
integral field spectroscopy with an instruments 
such as NIFS (McGregor et al. 2003) or SINFONI 
(Eisenhauer et al. 2003). 
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3.4. The low-mass end 

Out negatively curving Tlfbh " relation suggests 
that the SMBHs in dwarf elliptical galaxies will be 
smaller than predicted from the linear Mbh-»^ re- 
lation. When n = 1 {Mb ~ -14.5 ±1.5 mag; 
Graham & Guzman 2003, their figure 10), one 
can expect a SMBH mass 

and SAttl x 10^ using the quadratic and lin- 
ear relations, respectively. At n = 0.5 one would 
expect a SMBH mass of x IO^Mq and 

5.2^3'^ X IO^Mq from the quadratic and linear 
relations, respectively. When n = 0.25, the (ex- 
trapolated) quadratic relation predicts a mass of 
only ~ 2Mq. The absence of lO^M© black hole 
detections in local dwarf galaxies (e.g., Valluri ct 
al. 2005) and M33 (Gebhardt et al. 2001; Merritt, 
Ferrarese, & Joseph 2001) would argue against the 
extrapolation of the linear relation. The suitabil- 
ity of the curved log-quadratic Mbh-n relation to 
globular clusters (e.g., Gebhardt, Rich, & Ho 2005; 
De Rijcke, Buyle, & Dejonghe 2006, and references 
therein) remains to be explored. 

4. Central star clusters 

From an analysis of HST images of low- 
luminosity elliptical galaxies in the Coma clus- 
ter, Graham & Guzman (2003) presented a cor- 
relation between the luminosities of nuclear star 
clusters and that of their host galaxy. Similarly, 
strong correlations involving Virgo cluster ellipti- 
cals have been reported (e.g. Cote et al. 2006) and 
also between the nuclear star clusters in the bulges 
of spiral galaxies and their host bulge's flux (e.g., 
Balcells et al. 2003). The cluster-to- host flux ratio 
reported in Graham & Guzman (2003) was a few 
tenths of one per cent, intriguingly comparable 
to the mass ratio observed between SMBHs and 
their host bulge. In Figure 3 we have expanded 
the Mbh-Ti diagram by including the masses of 
nuclear star clusters plotted against their host 
galaxy's Sersic index. We have used two data sets 
obtained with the Hubble Space Telescope, and 
therefore likely to have had a reliable bulge/star- 
cluster decomposition. 

The first data set is from Graham & Guzman 
(2003, their Table 2). We have converted their 
F606W nuclear star cluster magnitudes into 
masses using an F606W-V color of —0.3 mag 
(Fukugita, Shimasaku & Ichikawa 1995) and a 



F-band mass-to-light ratio M/L = 1.5(M0/L0). 
This ratio is the median of the (King model) dy- 
namical M/L ratios from 57 star clusters listed 
in McLaughlin & van der Marcl (2005, their Ta- 
ble 13; see also Pryor & Meylan 1993, their Ta- 
ble 2). The second data set comes from Fer- 
rarese et al. (2006a, their Table 3), for which 
we have used an F475W-V color of 0.2 mag and 
M/Lv = 1.5{Mq/Lq). Due to the association 
of SMBHs with their host bulge, rather than the 
host galaxy, we have exluded the lenticular galax- 
ies^ because Ferrarese et al. (2006a) modeled the 
combined bulge+disc light with a single Sersic 
function. That is, no bulge/disc decomposition 
was performed. 

Figure 3 reveals (albeit with some scatter which 
is in part likely due to our conversion from fiux 
to mass) that the nuclear star clusters appear to 
follow the same log-quadratic relation as defined 
by the (Mbh,«) data set. Clearly there is no 
sharp transition; both SMBHs and star clusters 
have overlapping masses in the range 3 x 10^ < 
Miaco/MQ < 10^. Furthermore, the lower limit 
to the SMBH masses may be due to selection ef- 
fects, i.e., a reflection of our inability to resolve 
the SMBH sphere of influence in low-mass bulges 
(Merritt & Ferrarese 2001a). Similar results us- 
ing host bulge mass, rather than Sersic index, 
have recently been reported in Cote et al. (2006), 
Ferrarese et al. (2006b), Rossa et al. (2006) and 
Wehner & Harris (2006). McLaughlin, King, & 
Nayakshin (2006) propose that feedback from su- 
pernovae and stellar winds from a nuclear star- 
cluster can regulate the growth of the host bulge 
and thereby explain the latter connection. 

The result shown in Figure 3 is intriguing be- 
cause it suggests that the formation mechanisms 
of nuclear star clusters may share some common 
processes with the formation of SMBHs. The shal- 
lower potential wells, lower central stellar densi- 
ties (prior to core-depletion in massive elliptical 
galxies) and shallower central gradients in lower n 
bulges (Terzic & Graham 2005) may result in the 
production of star clusters rather than SMBHs. 

An obvious question is: Do any of the galaxies 
in Table 1 possess nuclear star clusters in addition 
to their SMBHs? For the Milky Way the answer 



Inclusion of the lenticular galaxies increases the scatter but 
does not suggest any shift in one particular direction. 
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is yes (Mcrritt 2006 and references therein). For 
the local group galaxy NGC 221 (M32), the excess 
central fliax above the Sersic model may not rep- 
resent a distinct additional central massive object, 
as is the case with the nuclei in NGC 205 (Valluri 
et al. 2005). In many galaxies, the presence of non- 
thermal flnx from AGN can complicate the identi- 
fication of nuclear star clusters when using only 
images and light-profiles. One additionally re- 
quires spectral information. In general, corc-Sersic 
galaxies do not show any excess nuclear flux, and 
the few which do show excess flux have an AGN. 
In contrast, the majority of the lower-luminosity 
Sersic ( "power- law" ) galaxies in Table 1 do show 
excess nuclear flux. However, this is usually due 
to an AGN. Two exceptions arc NGC 3384 and 
NGC 7457, which possess both a SMBH and a 
central star cluster (Ravindranath et al. 2001). 

While the mass of the star cluster in NGC 3384 
is approximately twice its SMBH mass, the star 
cluster in NGC 7457 is an order of magnitude 
greater than its SMBH mass* Creating a new 
Afbh-i-nico-'^ diagram would move NGC 7457 
(which has the second lowest SMBH mass in our 
sample) so that it overlaps with a cluster of three 
other data points seen in Figure 2 — leaving it 
consistent with the quadratic shown there. 

Successive mergers of nucleated, low-n, low- 
luminosity galaxies will eventually result in a 
galaxy whose mass is sufficient to expect a SMBH. 
This then raises the tantalising possibility of the 
recent formation of some SMBHs from the merger 
of two, or more, nucleated galaxies and the sub- 
sequent merger and 'runaway collision' of their 
respective nuclear star clusters. To date, such a 
coUisional process in star clusters has only been 
invoked to explain the formation of intermediate- 
mass black holcis (IMBHs: Qiiinlan & Shapiro 
1990; Portcgics Zwart et al. 1999; Freitag, Giirkan, 
& Rasio 2006). We are, however, unaware of the 
actual dwarf-dwarf galaxy collision rate today, 
and do not wish to suggest it is a frequent phe- 
nomenon. 

The formation scenario proposed here for po- 
tential young SMBHs is perhaps, at face value, 

'The apparent ff-band magnitudes were converted into ab- 
solute magnitudes using the distances in Table 1, and then 
into units of solar luminosity using MQ_jy=3.32 (Bessell, 
Castelli, & Plez 1998, their Appendix C and D), and fi- 
nally solar mass using M/Lh = ^■0{Mq/Lq). 



contrary to the conventional scenario in which 
they form via gas accretion during the height of 
quasar activity. Indeed, massive black holes have 
long been known to exist at high-redshifts, z > 3. 
Nonetheless, we do not rule out a kind of "down- 
sizing" or late-time formation of low-mass SMBHs 
built through a relatively dry runaway merger of 
nuclear star clusters. While involving no signif- 
icant AGN feedback, such a process could still 
maintain the Mmco/Afbuige mass ratio. 

// there are appreciable numbers of young 
SMBHs built from the merger of nuclei, this may 
require a downward revision to measurements of 
the average radiative efficiency of SMBHs and thus 
their implied rate of rotation (e.g., Ferrarese 2002; 
Shankar et al. 2004; Wang et al. 2006). 

5. Dependence on Hq 

The above linear and curved relations are 
nearly, but not completely, independent of the 
Hubble constant. While the Sersic indices arc in- 
dependent of galaxy distance, the SMBH masses 
depend linearly on the distance to each galaxy. 
Distances for 24 of our 27 galaxies were obtained 
from surface brightness fluctuation measurements 
and are thus independent of the Hubble constant, 
as of course is the Milky Way center. However, a 
Hubble constant of 73 km s~^ Mpc^^ was used to 
derive the distances to NGC 6251 and NGC 7052 
from their redshift. 

Studies which use the Mbh-n. relation, or in- 
deed the Mbh-cr relation, need to know if the 
SMBH masses predicted from the relation are de- 
pendent on Hq. If they are, then the Hubble 
constant needs to factored into such mass deriva- 
tions and thus also estimates of the SMBH mass 
function and space density. We have therefore 
rederived the quadratic Mbh-Ti relation excluding 
NGC 6251 and NGC 7052, so as to remove any 
possible dependency on Hq. Doing so, we obtain 
log(Mbh) = 7.97(±0.09) -I- 3.75(±0.46) log(n/3) - 
2.92(±0.83)[log(n/3)]2, with e = 0.17to o6 dex. 

This is not signiflcantly different from the ex- 
pression given in equation 7. In fact, the main 
reason for the slight shift in numbers is because 
of the exclusion of NGC 6251 — which resides at 
the high-n end of the relation, and thus has more 
weight than individual points in the middle of the 
relation — rather than the use of Hq = 73 km 
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s ^ Mpc ^ for this galaxy. This is evidenced by 
using Ho = 100 km s"! Mpc"^ for NGC 6251 
and NGC 7052 and adjusting their SMBH masses 
accordingly, which results in log(Mbh) = 7.97 + 
3.641og(n/3) — 3.11[log(n/3)]^. This expression is 
closer to equation 7 than the expression in the pre- 
vious paragraph is to equation 7, and reveals that 
the exclusion of NGC 6251 and NGC 7052 intro- 
duces more of a change than the use of an uncer- 
tain Hubble constant. It therefore makes sense to 
use equation 7 and, given the uncertainties on the 
fitted parameters a, b and c, treat equation 7 as if 
it were independent of the Hubble constant. The 
same statement is true for equations 2, 4 and 5. 

6. An Mbh— Mo relation 

There are now several relations between SMBH 
mass and the properties of the host bulge, such as 
luminosity, mass, velocity dispersion within some 
radius, and concentration. Which of these is the 
most fundamental remains unanswered (Novak et 
al. 2006). It may be that none of these are the 
drivers of the SMBH-bulge connection, but that all 
are symptomatic of a more profound relation. For 
instance, the combination of concentration and 
central density might be what is important — as 
this defines both the baryonic distribution and the 
strength of the gravitational potential with radius, 
at least until dark matter becomes a significant 
factor. It is therefore not unreasonable to expect 
that this might influence, or possibly even dictate, 
the ability of a bulge to fuel any central SMBH. 

The motion of the stars within a bulge, traced 
through the observable a, is the dynamical re- 
sponse to the underlying mass distribution. While 
a is therefore a tracer of mass, it is obviously de- 
pendent on how that mass is distributed, and is 
thus a) a function of the central stellar density 
and b) has a radial dependence which is set by 
the stellar concentration, at least within ~ li?e 
(Lintott, Ferreras, & Lahav 2006). 

While the SMBH mass might be a product of 
the radial mass distribution, specified by the cen- 
tral density and concentration, the observed re- 
lation between central density and concentration 
(e.g., Graham & Guzman 2003, their Figure 9; 
Merritt 2006, his Figure 5) subsequently means 
that one can express the SMBH mass in terms 
of just one of these quantities. Extrapolating 



under the nuclear star clusters in low-luminosity 
bulges, and over the partially-depleted cores in gi- 
ant galaxies — whose central stellar density has 
been modified by the slingshot effect of coalescing 
SMBHs — gives an estimate of a bulge's (original) 
central stellar density. Given the Mbh-ri relations 
derived in this paper, we can predict a correlation 
between Mbh and central stellar density. 

Substituting log(n) = (22.8 — ^o,b)/14, where 
fio,B is the central B-band surface density (Gra- 
ham & Guzman 2003)^, into equation 2 gives 

log(Mbh) = 10.91 - 0.19mo,b. (9) 

A relatively large uncertainty on the central sur- 
face brightness of I mag arcsec"^ translates to an 
uncertainty of only 0.19 dex in the logarithm of 
the SMBH mass, a value comparable to the in- 
trinsic scatter found for equation 7. Substitution 
of the above term for log(n) into equation 7 gives 

log(Mbh) = 8.13 + 0.25mo,b - 0.016(mo,b)^- (10) 

These two expressions differ most in their pre- 
diction of SMBH masses at the low mass end 
(Figure 4), with the linear relation predicting 
~ IO^Mq masses if /io = 25i3-mag, and the 
quadratic relation predicting masses of ~ lO'^M© 
at this faint central surface brightness. 

Although we are unable to conclude that an 
Mbh- Ho relation is more, or even equally, funda- 
mental than any other relation involving SMBH 
mass, we do consider it prudent to explore this 
possibility. We hope to acquire calibrated, high- 
resolution, near-infrared images to further inves- 
tigate this proposed relation. Due to the random 
projection angles of triaxial bulges on the plane of 
the sky, one might expect a certain amount of vari- 
ability in the observed values of fJ,o,B- We there- 
fore note that it may be advantageous to construct 
a plot of SMBH versus the deprojected (internal) 
central density. Such a program, however, is be- 
yond the intended scope of this paper. 

7. Summary 

We have constructed the Mbh-'fi relation us- 
ing updated SMBH masses and refined galaxy 



Many of our images arc photometrically uncalibratcd, so 
we are unable to plot log(Mbh) versus /io with the present 
data set. 
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Fig. 3. — Central 'massive compact object' mass 
versus host galaxy Sersic index 'n'. Filled circles 
represent SMBHs in elliptical galaxies, while open 
circles represent SMBHs in disk galaxies. The 
stars represent nuclear star clusters in elliptical 
galaxies presented in Graham & Guzman (2003) 
and Ferrarese et al. (2006). The straight line is the 
same as that shown in Figure 1, and the curved 
line is the same as that shown in Figure 2. 
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Fig. 4. — Predicted correlation between SMBH 

mass and the central B-band surface brightness of 
the host bulge, after allowing for excess nuclear 
flux and partially-depleted cores, as mentioned in 
Section 6. The solid line corresponds to equa- 
tion 9, the dashed line corresponds to equation 10. 



light-profile shapes, n. We performed an ordi- 
nary least squares, linear regression of Mbh on 
n that allows for measurement errors and intrin- 
sic scatter e, to obtain log(Mbh) = 7.81 (±0.08) -|- 
2.69(±0.28)log(n/3), with e = O.Utom dex in 
log Mbh- 

We additionally considered the possibility that 
the Mbh-'i' relation may not be linear, and applied 
a quadratic fit to account for the apparent curva- 
ture in the distribution. Although an arbitrary 
empirical function, the second order term in the 
quadratic fit 

log(Mbh) = 7.98(±0.09) + 3.70(±0.46)log(n/3) 
-3.10(±0.84)[log(n/3)]^ 

with an intrinsic scatter e = 0.18lQ;Qg dex, is in- 
consistent with a value of zero at the 99.99% confi- 
dence level. These relations were constructed with 
the objective of obtaining a relation that can be 
used to predict SMBH masses in other galaxies 
for which the bulge Sersic index is known. Un- 
like the linear relation, the quadratic relation pre- 
dicts finite SMBH masses at the high-n end, and 
^ mass black holes when n ^ 0.5, rather 

than ~ 1O^M0 obtained with the linear relation. 

The Sersic index, a quantity obtained from un- 
calibrated images, and independent of galaxy dis- 
tance, offers an easy way to acquire accurate es- 
timates of black hole masses in other galaxies. 
Moreover, the strength of the correlation (r = 
0.88, fs = 0.95) implies a fundamental connection 
that theory is yet to explain. In addition, we have 
shown that the nuclear star clusters at the centers 
of low- luminosity bulges appear to follow the same 
log-quadratic M^^h-n relation. This is suggestive, 
although not conclusive, that a similar formation 
mechanism may be responsible for the build up of 
massive compact objects (whether black holes or 
star clusters) at the centers of bulges. 

Finally, in Section 6 we have derived, for the 
first time, the expected linear and quadratic re- 
lation connecting the mass of a SMBH with the 
central surface brightness of its host bulge (Fig- 
ure 4). 

We arc grateful to Rebecca Koopmann who 
kindly supplied us with data for NGC 4649 (Koop- 
mann, Kenney, & Young 2001), and to Helmut 
Jerjen who provided light-profiles for NGC 3115, 
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profile for NGC 4564 was originally obtained by 
Peter Erwin as a part of Graham et al. (2001). 
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the statistical advice they provided early on in 
this work. We also thank Stuart Wyithe for his 
comments and kindly supplying us with a code 
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8. APPENDIX A: Galaxy light-profiles 

Wc have made \iso of the i?-band galaxy im- 
ages of NGC 4486, NGC 4564, NGC 4649, and 
NGC 4697 taken by Cheng et aL (1997), Frei et 
al. (1996), Kuchinski et aL (2000), and HST/ACS 
archive Prposal ID 10003 (PI C.Sarazin), respec- 
tively. In the case of NGC 3115 we used an /-band 
image from Kuchinski et al. (2000), and for the 
Milky Way we used the near-infrared light-profile 
from Kent, Dame, & Fazio (1991). 

For the ii-band images, the light-profiles were 
obtained in the standard way. Foreground stars 
and background galaxies were carefully masked, 
and the sky-background flux was determined from 
the mean of '--^10 median fluxes that were obtained 
from small boxes we positioned near the (galaxy- 
free) corners of each chip. The light-profiles were 
extracted using the IRAF task ELLIPSE, with the 
isophotal position angle and ellipticity free to vary, 
but the centroid fixed. 

An interesting deviation from the above recipe 
pertains to the analysis of the edge-on SO galaxy 
NGC 3115. Its light-profile was obtained us- 
ing another isophotal fitting routine written in 
IRAF (Jerjen, Kalnajs & Bingelli 2000). After 
the foreground and background objects, and non- 
symmetrical features about the nominal galaxy 
centre, were removed from the image, a symmet- 
rical two-dimensional (2D) model was constructed 
from the remaining light distribution by allowing 
the isophotal ellipticity and position angle to vary, 
but keeping the luminosity-weighted centre fixed. 
This process was repeated iteratively until the 
residuals were minimised. This left a residual im- 
age displaying only the edge-on disc of NGC 3115. 
The one-dimensional surface brightness profile of 
the (symmetrical) bulge was then calculated from 
this 2D model. 

Foreground stars were fitted with a Moffat func- 
tion that was used to quantify the point spread 
function in each image. However, in practice this 
wasn't important because in avoiding the poten- 
tial presence of nuclear star clusters or partially- 
depleted cores, we first removed the inner ~2 see- 
ing disks from the galaxy light-profiles before fit- 
ting Sersic's (1963) function. 

Each light-profile's best-fitting Sersic i?^/" 
function plus, in the case of NGC 4564, an ex- 
ponential disk, was obtained using the subroutine 



UNCMND from Kahaner, Molcr & Nash (1989). 
At each iteration, the nonlinear (seeing-convolved) 
Sersic function (plus exponential function when 
modelling NGC 4564) was approximated by a 
quadratic function derived from a Taylor series. 
The quadratic function was minimised to obtain a 
search direction, and an approximate minimum of 
the nonlinear function along the search direction 
found using a line search. The algorithm computes 
an approximation to the second derivative matrix 
of the nonlinear function using quasi-Newton tech- 
niques. The galaxy light-profiles and best-fitting 
models arc displayed in Figures 5 and 6. One can 
see, particularly from the residual profiles in the 
lower panels, that the fitted models perform well 
at matching the curvature in the light-profiles. 
This curvature, specifically, the Sersic index n, 
is observed to correlate strongly with a bulge's 
central supermassive black hole (e.g.. Figure 2). 
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Fig. 5. — Major-axis i?-band light-profiles for three elliptical galaxies with directly measured SMBH masses. 
The best-fitting Sersic function is shown by the solid curve. The arrows mark the spheroids' half-light radii. 
The residuals of the data about the fit are shown in the lower panels. (Note: Due to the fact that the Mbh-« 
relation does not require calibrated images, the profiles have not had their photometric zero-point determined.) 
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Fig. 6. — Major-axis light-profile for the disk galaxy NGC 4564, the bulge of the edge-on lenticular galaxy 
NGC 3115, and the bulge of the Milky Way. The best-fitting Sersic function (and exponential) are shown by 
the curved lines (and the straight lines). The solid lines show the seeing-convolvcd functions; the dashed lines 
show the unconvolved functions, i.e., the intrinsic profile. (Note: Due to the fact that the Mbh-n relation does 
not require calibrated images, the profiles have not had their photometric zero-point determined.) 
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